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The equivalence between thermal field theory and the Boltzmann transport equation is investi¬ 
gated at higher orders in the context of quantum electrodynamics. We compare the contributions 
obtained from the collisionless transport equation with the high temperature limit of the one-loop 
thermal Green’s function. Our approach employs the representation of the thermal Green’s func¬ 
tions in terms of forward scattering amplitudes. The general structure of these amplitudes clearly 
indicates that the physics described by the leading high temperature limit of quantum electrody¬ 
namics can be obtained from the Boltzmann transport equation. We also present some explicit 
examples of this interesting equivalence. 
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I. INTRODUCTION 


This paper is about hard thermal loop (HTL) am¬ 
plitudes in finite temperature quantum electrodynamics 
(QED). This has been the subject of many investigations 
not only in QED but also in non-Abelian gauge theories 
and gravity. These HTL amplitudes have loop momenta 
of the order of the temperature, T, which is considered 
large compared with all the external momenta. In the 
case of non-Abelian gauge field theories, it is possible to 
sum all the amplitudes in terms of a closed form expres¬ 
sion for the effective action IDS!- Similarly, in gravity, 
closed form expressions for the static and the long wave- 
lenght limits of the effective action have been obtained 
[4j [5]. The main motivation for these investigations is 
the search for a consistent thermal held theory in the 
perturbative regime. 

It is well known that the thermal effective action in 
QED has a leading temperature behavior proportional to 
T 2 which arises only from the contribution of the photon 
polarization tensor. All the other thermal amplitudes 
are either subleading or they vanish as a consequence of 
the neutrality of the plasma. On the other hand, it is 
also known that the HTL amplitudes can be generated 
using the Boltzmann transport equation. This approach 
has been used in scalar, non-Abelian gauge theories and 

gravity EHB3- 

When we adopt the Boltzmann transport equation ap¬ 
proach, in the case of QED, we realize that there is 
ground for considering some more general possibilities for 
the thermal amplitudes (although restricted to the semi- 
classical limit). First, since the plasma does not need 
to be neutral, amplitudes with an odd number of ex¬ 
ternal photons can be generated. Second, amplitudes 
which are not simply proportional to a power of T (like 
the T 2 behavior of the two-photon amplitude) are also 
generated. (This is the key distinction between the non- 
Abelian and Abelian theories, since in the former, each 
amplitude generated by the Boltzmann transport equa¬ 
tion formalism is proportional to the same power of T.) 
Furthermore, all these contributions are generated from 
a manifestly gauge invariant effective action, simply be¬ 


cause, as we will see, the Boltzmann transport equation 
formalism yields an effective action which is expressed in 
terms of the electromagnetic field tensor F^ v . 

However, it is not known if the same n-photon am¬ 
plitudes can also be obtained from a first principle ap¬ 
proach, like QED at finite temperature. Our main goal 
in the present work is to investigate this issue using the 
imaginary time formalism PEU. The main motivation of 
this investigation is to extend the concept of HTL ampli¬ 
tudes including all the contributions which arise from the 
Boltzmann transport equation approach. In fact, consid¬ 
ering that the Boltzmann transport equation formalism 
is a classical limit, it is important to understand how it 
relates with the high temperature limit of QED at finite 
temperature. As a first step toward this goal, we inves¬ 
tigate in the present paper the n-photon amplitudes in 
both formalisms. 

In the next section, we review the Boltzmann trans¬ 
port equation approach and derive the effective action 
and the corresponding thermal amplitudes. We obtain 
an implicit expression for all the thermal amplitudes and 
compute some explicit examples. In Sec. Ill we consider 
the thermal Green’s functions in QED at finite tempera¬ 
ture. We review how the thermal Green’s functions can 
be expressed in terms of forward scattering amplitudes. 
Then, using the HTL limit, we argue that the leading 
behavior of each of these forward scattering amplitudes 
agrees with the results obtained in Sec. [TT] The agree¬ 
ment is explicitly verified up to the four photon function. 
Finally, in Sec. |IV| we discuss the main results. 


II. THERMAL AMPLITUDES FROM THE 
BOLTZMANN TRANSPORT EQUATION 

Let us consider the semiclassical description of a QED 
plasma in terms of a distribution 

f(x, k , e) = e(k 0 )S(k 2 - m 2 )F(x , k , e), (2.1) 

where x = (t,x) and k = ( ko,k ) are the four- 

vectors of position and momenta. The Dirac delta 
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and theta functions impose the conditions k = 
(■ mjyj\ — v 2 , mv/y/l — v 2 ) for on-shell, positive energy, 
particles of charge e and mass to. 

Let us now assume that the collisions between the par¬ 
ticles are neglectable. Then, when an external electro¬ 
magnetic field is present, the trajectory followed by the 
particles is determined by the equations of motion 


to—— = m—— = eF^k ", (2.2) 

dr dr 

where F^ v = — d„A M is the electromagnetic field 

tensor and r is the proper time. As a result, the distri¬ 
bution function satisfies the equation 


d 

dr 


f(x,k,e) 


dj_ dx^_ df_dk E _ 

cb+ dr dk^ dr 


(2.3) 


We remark that the right-hand side of the previous equa¬ 
tion does not vanish when the effect of collisions is taken 
into account. Using Eqs. (2.2), we obtain the Boltzmann 
transport equation 


k ■ df(x, k, e) = -F • d k f(x, k, e), (2.4) 


where d k is the partial derivative in relation to the four 
momentum k a and we have defined 

T a = eF af5 kP. (2.5) 


This covariant form of the Boltzmann transport equation 
has been previously obtained in more general cases of 
non-Abelian gauge fields and also noncommutative gauge 
theories umm. 

Let us now consider the current which is produced 
when the plasma interacts with the external field. This 
can be expressed in terms of the distribution functions for 
the negative or the positive charge carriers (e.g. positrons 
or electrons) as 


= ±eC 


d A k k^fix, k, ±e). 


( 2 . 6 ) 


where C takes into account the degrees of freedom asso¬ 
ciated with the charged particles (C = 2 for an electron 
in 3+1 dimensions). The total current in a plasma would 
be given by 


J 11 — JU + Ju- 


(2.7) 


Of course one could think of other possible scenarios such 
that the net charge of the plasma does not vanish. For 
this reason, it is convenient to consider ± J fl separately. 
If we now perform the replacement k —» — k in Eq. 


(2.4), we obtain 


This shows that we can choose f(x , k, —e) = f(x, —k, e). 
Therefore we can, alternatively, express Eq. (2.6) as 


+» = ±eC 


J d 4 kk^f(x,±k,e), (2.10) 


so that the charge conjugation is equivalent to reversing 
the sign of the thermal particle momentum. 

Let us now consider the general relation that must exist 
between the current and the thermal effective action of 
electromagnetic fields interacting with a thermal plasma. 
This relation can be simply written as 

^transp = J d*X %(A) A» (x), (2.11) 


where we have made explicit the dependence of ,/„ on 
the external field which follows from the substitution 
of the solution of Eq. (2.4) into Eq. (2.10). 

In the most trivial scenario, when the external held is 
switched off, and the system is in thermal equilibrium, 
Eq. (2.4) is satisfied by 


f(x,k,e) = f (0 Hk) = j^26(k 0 )S(k 2 - m 2 )F^(k 0 ). 

( 2 - 12 ) 

The equilibrium distribution F'°'(ko) depends on the 
statistics of the particles. In the case of electrons and 
positrons 


F (0 \k 0 ) = N F (k 0 ) 


1 

exp ( k 0 /T) + 1 


(2.13) 


which is the Fermi-Dirac distribution (for our present 
purposes F^ 0 \ko) could be any ^-independent distribu¬ 
tion). Then, the partial currents are given by 



eC j" d 3 k 
(2?r ) 3 J \J\k\ 2 + to 2 

N f f \J\k\ 2 + to 2 ^ k^ (2.14) 


and the corresponding effective actions are 


ibp(l) 

1 transp 


± 


(2t r) 3 




M"(*)- 


(2.15) 


Introducing the four-vector K = (y 1 + m 2 /\k\ 2 , k/\k\), 
we can write 


k ■ df(x, -k, e ) = +F ■ d k f(x, -k, e). (2.8) 

On the other hand, if we make e —> —e we get 

k ■ df{x, k, -e) = +F ■ d k f(x, k, -e). (2.9) 


ipU) _ 4- 
transp 


CeT 3 1 (1 \m) 

( 2^)3 

d A x J dnK^ A'-' ix), 


(2.16) 
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where f dfl denotes the integration over the directions 
of K = k/\k\ and we have introduced the integration 
operator 


/W(m) = 


,4—r 


du 


1 


\Ju 2 + TO 2 e v/,,2 +™ 2 + 1 ’ 


TO 

m = 


(2.17) 


which operates on functions of ivo = i/l + m 2 /u 2 . (to¬ 
gether with the powers of the temperature, (to) give 
the full temperature dependence for the contribution of 
order n). Since f dflKi = 0, we finally obtain 


±p(i) = + eCT 3 

transp ^ 2^2 



J d 4 xA°(x). 


(2.18) 


Of course, in the case of a neutral plasma, the sum of the 
two partial contributions vanishes 


■p(i) _+p (i) 

1 transp 1 transp 


-r (1) = o 

transp 


(2.19) 


As we will see this simple property will remain true for 
all odd powers of the external field (see Eq. (2.25) ). 


In the presence of a general external field, the solu¬ 
tion of the Boltzmann transport equation can be quite 
involved. Here we consider external field configurations 
which can be dealt with using a perturbative approach. 
In these cases a formal solution of the Eq. (2.41 can be 
written as 

f( n \x,±k,e)=U 1 ^T-d k ' S ) f«>\k), (2.20) 

so that each /(") is of order n in powers of the external 
field and the complete phase space distribution can be 

expressed as /(x, fc,e) = f^° \k) + f^(x, k,e) +_ It 

is convenient to write Eq. ( |2.20[ ) in a more explicitly 
iterative fashion as 

/<">(*,±fc,e) = T e 

( 2 - 21 ) 

where n = 1, 2, .... From these solutions, and using 
Eqs. ( |2. 10 ) and (2.11), we can obtain the corresponding 
partial effective actions, defined order by order in terms 
of the phase space distributions, as follows 


±p( ra ) _ ± e r 

transp 


’ J d 4 x J d 4 k k/j, f( n 4 \x, ±k, e) 

( 2 . 22 ) 

In principle, the effective action is completely determined 
by this direct procedure starting from the knowledge of 
which can be chosen, for instance, as given by Eqs. 
( (2~12| ) and pH] ). 

The nontrivial effects of the external held interacting 
with the plasma will be manifest when we consider the 
second or higher order contributions in Eq. (2.22) (n > 
2). To derive the explicit dependence on the external 


Held, let us look at the details of the integrand in Eq. 
(2.22). Substituting Eq. (2.21) into Eq. (2.22) we obtain 


±p( n ) _ 2 n 

-tr'nnart ^ 


transp 


d 4 kA^k^^ZF a0 


—/ ( "- 2) (x,±fc,e). 


(2.23) 


Performing an integration by parts in the momentum in¬ 
tegral, yields 


±p( n ) 

1 transp 


= e 2 C / d 4 x / d 4 kA» 5 1 


1 


k^d a 


»k-d ( k-d ) 2 

k p F aP f^ 

= e 2 C J d 4 x J d 4 kA^ (5«k x d x - k^d a ) 

j^k^F^- 2 ), (2.24) 

where the derivative of d/dk a acting on kP produces a 
vanishing contribution as a consequence of the antisym¬ 
metry of the electromagnetic tensor Held. Performing an 
integration by parts in the configuration space integral 
and using again the antisymmetry of the electromagnetic 
tensor field we obtain 


±p( n ) _ p2 n 

x transp ° ^ 


d 4 x / d 4 kF^ 


1 


{k ■ d) : 


l 


kf J , 1 k^ 2 ij Vll , 2 f^ } (x, ±fc, e) 


(2.25) 


(from this result one can easily verify that, for odd n , 

+p( n ) 'p( ra ) __ O') 

x transp ' x transp w /’ 

This form of the effective action encodes a number of 


interesting features. First, since Eq. ( 2.25 ) is expressed 
directly in terms of F^ v (notice that is implicitly 

dependent on F^ as a solution of (2.21)), gauge invari¬ 
ance is explicitly manifested at any given order. Second, 
it is straightforward to obtain any higher order contribu¬ 
tions by successive substitutions of Eq. (2.21) into Eq. 


(2.25). These properties can be concisely summarized by 


expressing the external fields in terms of Fourier compo¬ 
nents so that the effective actions can be written as 


ir-.fi 1 ) _ _ 

1 transp 


1 f d 4 pi d 4 p 2 d 4 p n 


(2t r) 4 (2 tt) 4 ■ ‘ ‘ (2t r) 4 

A^( Pl )...A^( Pn ) 

(2tt) 4 5(pi H-h Pn), 


(2.26) 


where (PhP- 2 , ■ ■ ■ ,Pn) are the thermal ampli¬ 

tudes associated with the Boltzmann transport equation 
formalism. Then, from the gauge invariance of 
under A + Apf ! , we conclude that the thermal 
amplitudes must satisfy the Ward identities 


(pi,p 2 ,... J p n ) = o. 


(2.27) 
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Notice that in the case of non-Abelian gauge 
theories the corresponding Ward identities 
would relate (pi,p 2 , ■ ■ ■ ,P„) with 

n M?^ P Mn-i (jPi,P 2 ,---,Pn-i) so that all the ampli¬ 
tudes would have the same leading high temperature 
behavior. On the other hand, in the Abelian theory, 

each contribution rlnsp to the effective action is 
individually gauge independent and the corresponding 
amplitudes are not all proportional to the same power 
of the temperature. This can be simply understood by 
power counting the momentum dependence of successive 
terms obtained from Eq. (2.21). Furthermore, from the 
general form of Eq. (2.25) one can easily see that all 


the amplitudes will be functions of degree zero in the 
external momenta p. 

As an example, let us now compute the explicit form of 
the two-photon amplitude (thermal self-energy). Making 
n = 2 in Eq. (2.25) one readily obtains 


= e 2 C d 4 x d 4 kF^ 


1 


■ transp 

kpi Il y 2 f 

Using Eqs. (2.12) and ( |2.13[ ) 


(k ■ d)' 


.pin 


(o) 


(2.28) 


V 2 > e 2 CT 2 I^(fn) f 4 

4- 1/pn/n on /\ o I IX A 


• transp 


(2tt) 3 

J dQF^ K ^^2 11 ' 2 F™. (2-29) 


Expressing the external fields in terms of Fourier compo¬ 
nents we obtain 


±p(2) _ 

transp 


e 2 CT 2 I^ 2 \m) f d 4 p 

J (27 t ^ 

/ dn-- ,(2 ' 30) 


where 


F^ Vi = - p^rj^) Ap. (2.31) 


Inserting (2.31) into (2.30) and comparing the resulting 
expression with (2.26), we obtain the following result for 
the two-photon amplitude 


d zrrtransp 


2 e 2 CT 2 I^\m) 
(2tt) 3 



Fjj [t:y -(- J\ ,/Pfl 
K ■ p 


J dn 

p 2 k^k u \ 

( K-p) 2 ) 


(2.32) 


This particular example of the two-photon function il¬ 
lustrates the general structure of all ro-photon functions, 
which, as we will see can be expressed in terms of an 
angular integral over the directions of K = k/\k\, the in¬ 
tegrand being a rank n tensor which depends on AT M as 
well as the external momenta. The same structure has 
been previously obtained for the leading high temper¬ 
ature limit of non-Abelian as well as noncommutative 
gauge theories m- 


Iterating Eq. (2.21) one more time, Eq. (2.25) yields 


±r'( n ) 

1 transp 


T e 3 C J d A x J d A kF^ k ^ 2 ^ 2 F^ 

( k ^ F ° e i^)’ in - 3 F±k,e). (2.33) 


Performing an integration by parts in the momentum in¬ 
tegral, the result can be cast in the form 


±p(«) _ ± e 3 r 

transp ° 


d 4 x / d 4 kF^ 


d Xl 


(.k ■ d) : 


2 


d X2 d X3 

(k ■ d ) 3 


F^ 3 3r h'n'2T l j, lll2fl3l/3 \ 1 \ 2 \ 3 


/ (n “ 3) 0 v,±k,e), (2.34) 


where 


r Fi_i 1 ^ 2 fi 3 u 3 \ 1 x 2 x 3 — ‘^kf 2l kf l2 kf l3 kx 2 k\ 3 if]x 1 v 3 [kmk/ J , 3 kx 1 kx 2 kx 3 1 )^ 2 ^ U Mi ^ M 2 ] + k lil kf J , 2 k /J , 3 kx 1 kx 2 ^ 3^3 (2.35) 


It is then straightforward to obtain the three-point am¬ 
plitudes. Similarl y to t he previous cases, we now consider 
n = 3 in the Eq. (2.34) and express the external fields in 


terms of Fourier components. Comparing the resulting 
expression with Eq. (2.26) we obtain the result 


drrr transp _i_ 

LL p 1 p 2 p3 ~ 


6e 3 CTI^(fh) 

(>rj3 


dn A MlM2M3 .(2.36) 


where the explicit result for A^ lfl2lJ , 3 is given in Appendix 

El 

Using this iterative approach, we have also computed 
the explicit result for the four-photon amplitude. Since 
the calculation is straightforward, but rather involved, 
we have made use of computer algebra. The result has 
been compared with the one obtained using the high tem¬ 
perature limit of thermal field theory, as we describe in 
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the next section. 

It is clear that, in general, the iterative procedure will 
produce a result for the ?r-photon function which is pro¬ 
portional to 

< 2 - 37 > 

so that each thermal amplitude has a simple temperature 
dependence given by the factors T A ~ n I^ n ' ) {fh). 

III. HARD THERMAL LOOPS IN QED 

In this section we consider the QED ra-photon ampli¬ 
tudes in the usual thermal field theory setting mm - 
Furthermore, we chose the imaginary time formalism, 
which, as it will be seen, is the most direct approach 
to obtain the thermal amplitudes in the present context. 
Also, to fully describe a plasma under the influence of 


an external electromagnetic fields, the external energies 
in these diagrams have to be analytically continued to 
continuous values, after the Matsubara sums have been 
performed. This prescription extends the imaginary time 
formalism to a nonequilibrium regime under the influence 
of external fields. 

Let us consider the one-loop diagram, shown in Fig. |T| 
Once we compute this basic diagram, the n-point one- 
particle-irreducible (1PI) amplitudes can be obtained as 
the sum of certain permutations of the pairs (pi, Pi), so 
that the bosonic symmetry is satisfied. In the imaginary 
time formalism the contribution of this basic diagram can 
be written as 

..^n = ~T f E iM2...^( fc o,fc) (3.1) 

J ' k 0 =ioj n 

where 


Xko,k) = 


S k ',Pl,-- ,Pn) 


k 0 - \k\ 2 - rn 2 (k 0 + Pi 0 ) 2 - (k +pi) 2 - m 2 (k 0 - p n0 ) 2 ~ (k - p n ) 2 ~ m 2 


(3.2) 


and we have taken into account the minus sign associated 
with the fermion loop. The factor (fc;pi, •• • ,p n ) 

is a shorthand notation for all the contributions to the 
numerator which arises, for instance, from the traces of 
Dirac matrices and coupling constant factors. The quan¬ 
tities Lo n = (2 n + 1)ttT are the fermionic Matsubara fre¬ 
quencies. Using the identity mm 


T E f(ko = 


iu n ) = dk 0 f(k 0 ) 


tanh ( ^/3 k 0 


(3.3) 


it is straightforward to show that the temperature de¬ 
pendent part of (3.1) can be written as 


G 


therm 

V>lP2—Pn 


d 3 k 


pioo+S 


S^(fco) 


’ ( 27r ) 3 J-ioo+S 2™ 

fII iAt2 -A»n (ko, k) + ( — ^Oi k) 


Making the change of variables k —> — k in the second 
term inside the bracket and using the shorthand notation 

fnui 2 ...iin(koi k) = finuz-.-nn {k)i we can write the thermal 
contribution of the n-photon diagram in terms of two 
components, as 


xn itherm _ -\-/^itherm 

L Vi/L 2 ...Ai Tl — ^ 


—s'itherm 


+ ~G 


••■Mn ’ 


(3.5) 


where 


±Qtherm 


f d 3 k r io ° +6 dk 0 

/ TTAs / ttwIVf fc 0 
J (27r) d J„ ioo+s 2 ?ti 

/ MlM2 ... M „(±fc). (3.6) 


This separation in two components is the first necessary 
step in order to compare the field theory formalism with 
the results obtained in the previous section. 


Closing the contour of the k 0 integration in Eq. (3.6) 
in the right-hand side of the complex plane and taking 


into account the poles from the denominator of Eq. (|3.2|) 
we obtain 


±£itherm 


d 3 k 


n f v ( fc ) 2 + ; 




(2t r) 3 


2G(k) 2 


m 


cycl. 


(3.4) where the quantities AJJ 


denote on-shell 

(ko = \J (k) 2 + m 2 ) forward scattering amplitudes, 
as depicted in Fig. [2j and J2 C yci denotes the sum of 
all cyclic permutations of the pairs ( pi , pi), as we will 
explain soon. 

At this point, we analytically continue the result 
to the Minkowski space-time, so that the amplitudes 
A" (k;pi,p2, ■ ■ ■ , Pn) can be computed using the 
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FIG. 1: One-loop amplitude which contribute to the effec¬ 
tive action. The solid and wavy lines represent fermions and 
photons respectively. The external momenta are denoted by 
Pi, i = 1,... n and k denotes the loop momenta. The sym¬ 
metrized amplitude which contributes to the effective action 
is obtained by adding all the permutations of (n— 1) external 
photon lines, and dividing the result by (n — 1)! . 


zero temperature QED Feynman rules. In this way, we 
obtain from Fig. [2] the following result (an extra minus 
sign arises from the analytic continuation) 


A l ± = -etr(7a 1 7Mi)fc“ 1 = ~4e/c Ml (3.8a) 
= _e " ( tr 7«i7Ml ■■•7a„7/uJ 
k ai ... (k + s n - i) a ™ 


(pf + 2k-pi)... (s 2 _i + 2k ■ s„- 1 ) 


, (3.8b) 


where 


s n -i = Pi + P 2 - hpn- 1 ; n > 2. 


(3.9) 


We remark that, as a consequence of the on-shell con¬ 
dition, all the even powers of m cancel. Also the terms 
proportional to odd powers of m vanish simply because 
they are proportional to the trace of an odd number of 
gamma matrices, which is zero when the space-time di¬ 
mension is even (in an odd-dimensional space-time there 
would be induced Chern-Simons terms proportional to 
the odd powers of rn). This way of expressing the 71- 
photon amplitudes in terms of traces of 2 n gamma ma¬ 
trices is convenient in order to carry a computer algebra 
calculation. 

The partic ular set of forward scattering amplitudes 
given by Eq. (3.8) includes the ones which arise from the 


residue of the pole of 1/fc 2 in Eq. (3.2). It is straightfor¬ 
ward to show that the corresponding contributions asso¬ 
ciated with each pole of Eq. (13.21) produce a result which 


is given by a cyclic permutation of the pairs ( pt , /i,;). In 
this way, the final result can be expressed as the sum 
E cyc z i n Eq. (3.7). This representation of one-loop 1PI 
diagrams in terms of forward scattering amplitudes of on- 
shell thermal particles has been derived in various other 
physical systems such as non-Abelian gauge theories and 
gravity [l8l - [21] . As a result, in all these cases it is nat¬ 
ural to think in terms of a physical scenario where on- 
shell thermal particles of momentum k are scattered by 


the external photons of momenta Pi, i = 1,..., n, as ad¬ 
vanced by Barton [22]. This is clearly in tune with the 
Boltzmann transport equation approach of the previous 
section. 

The full contribution to the 1PI thermal Green’s func¬ 
tions can now be obtained adding all the (n— 1)! permu¬ 
tations of the basic result in Eq. (3.7). Then, taking into 

= fZ 


account that E^™- E cy ci. 
the 1PI amplitudes as 


we can express 


rbrr therm 


d 3 k 


n f 


(ky 


J (2tt) 3 

E A , 


2 Y (fc) 2 + TO 2 

,Pn), 


(3.10) 


and the corresponding effective actions are 

±p(n) _ 1 f d 4 pi d 4 p 2 d 4 p n 

therm ~ (n — 1)! J (2tt) 4 (2tt) 4 " ' (2tt) 4 

A^(p 1 )A^{p 2 )...A^{p n ) 

± n^ i e ; 2 m .^(Pl,P2 ,...,Pn) 

(2tt) 4 <5(pi -I-b Pn). (3.11) 


(We remark that, due to the symmetry of the basic di¬ 
agram in Fig. |T] the sum of the (n — 1)! permutations 
of external legs produces completely bosonic symmetric 
amplitudes ± n^ i e ™ ( Pl ,... ,p n ).) 

In the pre vious section we have seen that all the effec¬ 
tive actions (2.26) have a simple temperature dependence 
of the form T 4- ”/( n )( to), as a result of Eq. (2.37). On 


the other hand, the thermal field theory effective actions 


(3.11) have a more complicated temperature dependence. 


From a physical point of view, we expect that the Boltz¬ 
mann transport equation approach is restricted to the 
high temperature limit, when the system approaches a 
classical limit. This indicates that in the thermal field 
theory description, we have to consider the limit when 
T is high compared with the energy-momentum scale 
of the external field. This amounts to consider the so- 
called hard thermal loop approximation, which in turn 
implies that the integration in (3.10) is dominated by 
the region where k^ ~^>Pi^. Denoting the HTL limit of 

±n ^i e ™>i< • • • b y ±n w...M> 1 > • ■ • equiv¬ 

alence of the two effective actions would then imply that 


±jjhtl 




> 


l 


ii •• • >j 


ft>) = - ±I C7.(P 1 .‘-.ft,). (3-12) 


In QED this is a well-known result for the leading T 2 
contribution to the two-photon amplitude. (In the case 
of non-Abelian gauge theories it can be proved that all 
the HTL n-gluon functions, which have the same T 2 be¬ 
havior, can be generated using the Boltzmann transport 
equation approach [23].) Our main goal here is to inves¬ 
tigate the validity of this equivalence also for all other 
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P 1 M-; P 2 M-2 P 3 M-i 



Pn F„ 


In the hard thermal loop approximation the denomina¬ 


tors of the forward scattering amplitudes in (3.13) can 
be expanded using 


k-P, 


FIG. 2: Diagram representing the forward scattering ampli¬ 
tude A^ 1M2 ... Mii (fc,pi, • • ■ ,p n ). A trace over the product of 
gamma matrices is to be understood. 


n-photon amplitudes, which have each a distinct depen¬ 
dence on T, in the HTL approximation. 

Using the explicit expression for the on-shell momen¬ 


tum = |fc|Ar M , with K = (yl + m 2 /\k\ 2 , k/\k\), Eq. 


(3.10) yields 


±1 I htl 


3 , \k\ 2 d\k\- 


N f (sj(k) 2 +m 2 


( 27r )- Jo 

r dn ^2 a l 


2 y ( k ) 2 + m 2 

piP2...n n ( ± -\k\ K ,Pi,-- ,Pn)■ 
(3.13) 


1 _ 1 1 _ 1 s 2 n 

s 2 + 2/c • s n |fc| 2K ■ s n |fc | 2 (2 1\ ■ s n ) 2 


+ .... (3.14) 


When combined with the momentum dependence of the 
numerators, as shown in Eq. (3.8), the resulting expres¬ 


sion for the forward scattering amplitudes in Eq. (3.13) 
can be expressed as 


Ml At 2 • • • Mi 


(±\k\K,pi, ■ ■ ■ ,p n ) = ±\k\A\ 


A i lA i 2---A i r 






±-^Al 
\k\ ‘ 




\kf 


-a: 


M1M2---M7 


zb, 


(3.15) 


r 


Let us first consider the simplest case, namely the one- 
photon amplitude. In this case the forward scattering 
amplitude is exactly given by 


A 1 (±uTK) = TeuT tr#^ = =f4 euTK^ 


(3.16) 


where uT = |fc|. Using Eq. (3.10) we obtain the following 


result for the partial thermal one-photon amplitudes 

T rm = ± 2er ^)3 (m) / dnKf *' (3J7) 

Comparing this result with the left-hand side of Eq. 
(3.12) (for n = 1) , which can be easily obtained from 
Eq. § T 6 |, we can see that Eq. (|3.12|) is satisfied pro¬ 


vided that C = 2, which is the correct result for fermions 
in 3 -T1 dimensions. 

Next, let us consider the more interesting cases when 
n > 2 . As we have seen in the previous section, 
the Boltzmann transport equation approach produces 
(n > 2 )-photon amplitudes which are proportional to 
r 4 -"/(")(m). On the other hand, Eq. (3.15) seems to 
yield, for each thermal amplitude, when substituted in 
Eq. (3.13), several other terms proportional to higher 


powers of T. Therefore, a necessary condition to have 
the equivalence ( 3. 12 [ ) , is that contributions proportional 
to higher powers of T should vanish. For instance, the 
term proportional to in the two-photon func¬ 

tion, as well as the terms proportional to T^I^ijn) 
and T 2 I^ 2 \i fi) in the three-photon function, and so on, 
should vanish. 

All the contributions proportional to T 3 /W (m), which 
would contribute to (3.13), come from the first term in 
Eq. (3.15). From the HTL procedure above described, 
one can see that A^ fii ^ ^ must have the highest degree 
in K, being proportional to 


K K K 

lv cx.2 • • • 1Y ct n 

{K ■ Pl) [K ■ (p! + p 2 )] ■ ■ ■ [K ■ (pi +p 2 + Pn- 1)] ’ 

( 3 ; 18 ) 

After performing the sum over all the cyclic permutations 
and using momentum conservation, it is straightforward 
to show that the resulting expression vanishes. More for¬ 
mally, we have the following necessary condition for the 
equivalence high-temperature QED and the Boltzmann 


































transport equation formalism of the previous section. 


of Eq. (3.21). Indeed, if we contract it with p 3 we 




,p n ) = 0; n > 2. (3.19) 


cycl. 


We remark that this is a property of individual dia¬ 
grams which is true at the level of the forward scatter¬ 
ing amplitudes (one does not need to add the (n — 1)! 
permutations to cancel these so-called superleading T 3 
contributions). Of course for a neutral plasma the sum 
+ n u!u 2 ..m + would vanish trivially, since the 

terms which are proportional to any odd power of T are 
also odd in k. However, in the present context, it is 
important to check the behavior of the individual com¬ 
ponents (for both even and odd n), in order to properly 
compare with the Boltzmann transport equation results 
of the previous section. 


Since Eq. (3.19) already eliminates the T 3 contribution 


to the two-photon function, only the second term in Eq. 
(3.15) may contribute to leading order, producing a T 2 


contribution. A straightforward calculation yields (p 3 = 
p and p 2 = —p) 


A 2 


o 2 / K-iyPfi P KilKv 

= 2e - 


(- K-p ) 2 

(3.20) 

Inserting Eq. (3.20) into Eq. ( |3. 13 ), we obtain 


±i | htl _ _ 


J dfl - 


2 e 2 T 2 /( 2 )(m) 


(27r) 3 

K,jp v + K u p^ p 2 K tl K v 


K ■ p 


(k-p¥ 


■ (3.21) 


Comparing Eq. (3.21) with Eq. 
(3.12) is satisfied for C = 2 


whic 


(2.32) we can see that 


r is the correct result 


for fermions in 3 + 1 dimensions. 

Up to this point the thermal field theory calculations 
have been performed without using any computer algebra 
tool. If we go further and try to verify Eq. (3.12) for 
n = 3,4,... the calculation starts to become much more 
involved. Therefore, to proceed up to at least n = 4we 
have employed the FeynCalc computer algebra tool [23] , 

At the three-photon order, there would be in princi¬ 
ple contributions from A%, which would give rise to a 
contribution proportional to T 2 . This would invalidate 
the equivalence of the Boltzmann transport equation ap¬ 
proach with thermal field theory. However, a straightfor¬ 
ward calculation yields 


M 


2e 3 A^ ; 
K ■ p 3 
(pi) 2 K t 




K 


{K-p i) 2 


_ K^iPi^ + K^pi^ 

K ■ Pi 

- (Hi,Pi) O- (// 3 ,p 3 ), 

(3.22) 


which is manifestly antisymmetric. Therefore, the sum 
of its permutations vanish trivially. (It is remarkable 
how Eq. (3.22) is simply related to the tensor structure 


obtain a relation which resembles the Lorentz structure of 
the non-Abelian gauge theory, without the antisymmetric 
colour factor). 

We are then left with a leading contribution to the 
three-photon function which is of order T as in the case 
of the Boltzmann transport equation approach. With 
the help of computer algebra, the resulting expression 
for Aj can be computed in a straightforward way. 

After performing the permutations of the external photon 
lines, we obtain complete agreement with (3.12), when 
using Eqs. (2.36) and (Al). 


The next and more interesting n = 4-photon order can 
be dealt with similarly to the previous cases, although the 
calculation becomes much more involved. A straightfor¬ 
ward computer algebra calculation shows that, besides 
the cancellation of the contribution proportional to T 3 , 
the other contributions due to Ai „ ,, andil, ,, also 

vanish after we add the permutations. Finally, the re¬ 
maining result is also in agreement with ( 3. 12[ ) as in the 
previous cases. This shows that, up to the four-photon 
order, we have full agreement with the Boltzmann trans¬ 
port equation approach. 

The previous results show that equivalence between 
thermal field theory and the Boltzmann transport equa¬ 
tion approach is verified at the level of the forward scat¬ 
tering amplitudes, in the HTL limit, so that it is not 
necessary to explicitly perform the angular integrations 
in Eq. (3.13). 


IV. DISCUSSION 

In this work we have investigated the equivalence be¬ 
tween the high temperature limit of thermal field theory 
and the Boltzmann transport equation. Specifically, we 
have argued that both approaches produce the same n- 
photon thermal amplitudes. This is indicated by the gen¬ 
eral form derived for the amplitudes, as well as from the 
explicit calculations up to n = 4. Furthermore, we have 
considered the more general cases when the plasma may 
not be neutral, so that the amplitudes with odd values 
of n are nonvanishing. This equivalence implies that the 
quantities A” in Eq. (|3.15|) such that p < n should van¬ 


ish. Only the contribution with p = n can be interpreted 
as arising from the classical limit of finite temperature 
QED. On the other hand, the p > n contributions repre¬ 
sent quantum corrections which are not described by the 
Boltzmann transport equation. 

A similar equivalence is well known in the case of non- 
Abelian gauge theories as well as gravity and also in the 
particular case of the two-photon amplitude in QED. In 
the case of non-Abelian theories it has been useful in or¬ 
der to derive a closed form for the gauge invariant hard 
thermal loop effective action [23] . Our analysis shows 
that all the results obtained from the Boltzmann trans¬ 
port equation can be identified with a more general HTL 
limit of the thermal Green’s functions in QED. Conse- 
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quently, the gauge invariant effective action, rather than 
being proportional to a single power of the temperature, 
will have a more general dependence on T as can be seen 
from the form of the thermal amplitudes in Eq. pl7| . 
We remark that the electron mass is not being neglected 
so that the amplitudes with n > 4 are well defined func¬ 
tions of m/T. 

Once the equivalence with the Boltzmann transport 
equation is established, a full resumation of the thermal 
effects in QED, to all orders in the external field, would 
require a further investigation, taking into account the 
gauge invariant effective actions as given by Eq. (2.251. 
The resulting effective action would be relevant to de¬ 
scribe a low density plasma in an external field with fre¬ 


quencies which are small compared with the temperature. 
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Appendix A 


Here we present the result for the integrand of the 
three-photon amplitude (2.361: 


An 


E 

perm. 
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2 K fl3 K, 
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2 K M3 AT Mi p 3fl2 Pi-P2 K1<^ p lfl2 p 2 ■ p 3 
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ix 2 P 3 tll P2^ 3 _ I<^ 2 p 2 -P 3^ iM 3 _ K ix 1 K ■P2Pl 2 r] l x 2 iX3 _ J Vl P 2 fx 3 P^ix 2 
)K • p 3 ) 2 ( K-p 3 f K-p 3 (K-pi) 2 K-psK-pi 

K IX 1 K ■ P2 P3i2 2 Plp 3 K -P 2 P 1 - P3 V1X2IX3 , K -P2 Pl^ Vix 21x3 


(■ K-P3 ) K-pi 
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I< ■ p 3 
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(- K-P3f 


{K ■ p 3 ) K Pl 

K -P 2 P 3 ^ V/xi M 3 

(■ K-P3f 


K ■ p 3 K ■ Pi 


r 


(Al) 


A similar but much more involved expression has been 
obtained for the four-photon amplitude, using computer 
algebra. In all the cases we have obtained agreement 


between the results derived using either the Boltzmann 
transport equation or the high temperature limit of ther¬ 
mal field theory. 
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